Levinson theorem for Dirac particles in one dimension* 

Qiong-gui Lin^ 

China Center of Advanced Science and Technology (World Laboratory), 
P.O.Box 8730, Beijing 100080, People's Republic of China 

and 

Department of Physics, Zhongshan University, Guangzhou 510275, 
People's Republic of China $ 



The scattering of Dirac particles by symmetric potentials in one di- 
mension is studied. A Levinson theorem is established. By this theorem, 
the number of bound states with even (odd) parity, n + (n_), is related 
to the phase shifts rj + (±E k ) [r]_(±E k )} of scattering states with the same 
parity at zero momentum as follows: 



rj ± (/i) +r]±(-fi) ± -[sin 2 ?7 ± (/i) - sin 2 r/± (-//)] = n±n. 



The theorem is verified by several simple examples. 



PACS number(s): 34.10. +x, 03.65.-w, 11.80.-m 

*published in Eur. Phys. J. D 7 (1999) 515-524. 
^E-mail addresses: qg_lin@163.net, stdp@zsu.edu.cn 
* Mailing address 

1 



I. INTRODUCTION 



In 1949, Levinson established a theorem in nonrelativistic quantum mechanics [1]. The 
theorem gives a relation between bound states and scattering states in a given angular mo- 
mentum channel /, i.e., the total number of bound states n\ is related to the phase shift Si(k) 
at threshold (k — 0): 

<J,(0)=n,7r, / = 0,1,2,.... (la) 
The case I = should be modified as 

5o(0) = (n + 1/2)tt (lb) 

when there exists a zero-energy resonance (a half-bound state) [2]. This is one of the most 
interesting and beautiful results in nonrelativistic quantum theory. The subject was then 
studied by many authors (some are listed in the Refs. [2-8]) and generalized to relativistic 
quantum mechanics [6,9-15]. In Ref. [6] the theorem was first written as 

*i(0) -(-)'! sin 2 5,(0) = ni 7r. (2) 

The second term on the left-hand side (lhs) appears automatically in the Green function 
approach to Levinson's theorem, which was used in earlier works [3,4] and further developed 
by Ni [6]. This method is quite different from the original one where the theorem is obtained 
by using the analytic properties of the Jost functions, and is convenient for generalization to 
the relativistic case. It can be shown that Si(0) /ir takes integers when / > or when I = but 
there exists no zero-energy solution. Then Eq. (2) is equivalent to Eq. (la). When / = and 
there exists a zero-energy solution (a half-bound state), however, 5i(0)/tc takes half integers 
and sin 2 <5 (0) = 1. In this case Eq. (2) is equivalent to Eq. (lb). Thus Eq. (2) combines 
the two cases in Eq. (1) into a unified form. The correct generalization of Levinson's theorem 
to Dirac particles was first obtained by the Green function method in Ref. [10]. If the phase 
shifts in the angular momentum channel k (k — ±1, ±2, . . .) are denoted by r) K (±Ek) and the 
total number of bound states in that channel by n K , then the Levinson theorem read 

vM+Vk(-v) ~ (-) K e(/t)-[sin 2 r] K (/i) - sin 2 rj K (-fj,)] = n K n, (3) 

where e(/c) = 1 (—1) for positive (negative) k, \x is the mass of the particle. The effect of half- 
bound states has been automatically included in this equation, though the problem is more 
complicated than in the nonrelativistic case. There is a modulo- 7r ambiguity in the definition 
of the phase shifts. This is common to both nonrelativistic and relativistic cases. In the former 
case this was resolved by setting £/(oo) = (rather than a multiple of n) which can be freely 
done. In the latter case one is not allowed to set r] K (±oo) = 0, but the modulo-7r ambiguity 
can be appropriately resolved [9, 10, 16 17]. 

Most of the above cited works mainly deal with the problem in ordinary three-dimensional 
space. Recently, several authors have studied the two-dimensional version of Levinson's theo- 
rem for both nonrelativistic Schrodinger particles and relativistic Dirac particles [18-20]. In the 
works of the present author the Green function method is employed [18]. In the nonrelativistic 
case we obtain 

5 m (0) = n m 7c, m = 0,1,2,.... (4) 



2 



In two dimensions S m (0)/n always take integers. This threshold behavior is quite different 
from that in three dimensions and in one dimension (see below). In Ref. [19], Dong, Hou, and 
Ma studied the problem by the method of Sturm-Liouville theorem and found that the case 
m — 1 in Eq. (4) should be modified as 

MO) = (ni + 1)tt (5) 

when there exists a half-bound state in this channel. In other words, the half-bound state 
with m = 1 plays the role of a real bound state in the Levinson theorem. The same result 
was also obtained in an earlier paper [20]. Similar modifications also occur in the relativistic 
case. The miscounting in our work may be due to the very subtle behavior of the above state 
and the shortcomings of the regularization procedure we employed [cf. Eqs. (44-47) in this 
paper], rather than the failure of the completeness relation as remarked in Ref. [19]. (see some 
more discussions in Sec. II.) The above state is quite different from the half-bound state with 
m = and from the half-bound states in three and one dimension because it tends to zero at 
infinity. As a result it is lost in the regularization procedure. In three- and one-dimensional 
space there is no similar case and the Green function method leads to correct results. 

It seems that the one-dimensional version of Levinson's theorem has attracted more atten- 
tion than the two-dimensional one. In fact, the nonrelativistic case has been studied by several 
authors [21-23]. In a symmetric potential V(x) (an even function of x), the Levinson theorem 
takes the form 

5 ± (0)±|sin 2 5 ± (0)=n ± 7r, (6) 

where + (— ) indicates even (odd) parity. It is easily seen that the odd-parity case coincides 
with the case I — in three dimensions, while the even-parity case has no counterpart in three 
dimensions. This is the main feature in one dimension. 

The purpose of this paper is to generalize the one-dimensional Levinson theorem to rela- 
tivistic Dirac particles. The interest of this problem is threefold. First, the Levinson theorem 
is a nonperturbative result. It may be useful in the study of nonperturbative field theories. In 
view of the wide interest in field theories and condensed matter physics in lower dimensions 
in recent years, and the potential applications of the Levinson theorem, the problem seems of 
interest. Indeed, some applications of the theorem to field theories have appeared in the liter- 
ature [21, 24, 25]. More recently, the Levinson theorem in two dimensions has been used [26] 
to the study of the screened Coulomb potential which plays an important role in the physics 
of semiconductor heterostructures. Second, as an academic problem, it exhibits some new 
features as compared with the three- or two-dimensional problem. Thus it may be of some in- 
terest in itself, as many one-dimensional models, say, the one-dimensional QED (the Schwinger 
model), studied in the literature. Third, it may be related to the real three-dimensional prob- 
lem. Consider a Dirac particle in an external potential V (the zero component of a vector 
potential). If V depends only on one space variable, say, x, one may consider a special case 
where the particle moves only in the x direction. This special case may be described by the 
one-dimensional Dirac equation. This connects the problem to the real physical world. 

Consider Dirac particles in a symmetric potential in one dimension. The phase shifts of 
even- (odd-) parity solutions (see Sec. II) are denoted by r)+(±Ek) [r)-(±Ek)], while the number 
of bound states with the same parity is denoted by n + The Levinson theorem connects 

them with each other: 
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V±(v>)+V±(-t*) ± -j[^ 2 ri±({i) - sin 2 ?? ± (-^)] = n±n. 



(7) 



We will establish this theorem for cutoff potentials (vanish when \x\ > a > 0) in this paper. 
Though the even-(odd-)parity case coincides in appearance with the case k — 1 (k — — 1) 
in three dimensions [cf. Eq. (3)]. They are essentially different. First, the one-dimensional 
equation is different from the radial one with k = f or k = -1 in three dimensions. The 
boundary conditions at the origin are also different. Second, the two critical energy solutions 
with E = ±/x and even (odd) parity, if exist, are both half-bound states in one dimension, 
while one of the two in three dimensions with k — 1 (k — — 1) is a bound state. Third, the 
threshold behaviors of the phase shifts are different. For example, in three dimensions rji(/j,)/ir 
always takes integers, while in one dimension r] + ( y fi)/iT takes integers only when there exists a 
half-bound state with E = ji and even parity. In general r) + ((j,)/ir takes half integers. There 
is a similar difference between rf-(—fi) and 77-1 (— fi). Therefore each case in Eq. (7) has no 
counterpart in three dimensions. 

Throughout this paper we use natural units where H — c — 1. In the next section we first 
discuss the solutions of the one-dimensional Dirac equation in an external symmetric potential. 
Then we discuss the scattering problem. In Sec. Ill the behavior of the phase shifts near k = 
is analyzed. The behavior when there exists a half-bound state is distinguished from the case 
without one. In Sec. IV we establish the Levinson theorem for cutoff potentials. In this paper 
we do not resort to the Green function method developed in Ref. [6] and used in Refs. [10, 17- 
18]. We employ the completeness of the solutions of the Dirac equation to derive the theorem 
more directly. This is essentially equivalent to the Green function method but simpler, and has 
been used in some recent works [13, 23]. In Sec. V, we first discuss the behavior of the phase 
shifts at infinite momentum and the resolution of the modulo- n ambiguity in the definition of 
the phase shifts. Then we examine the Levinson theorem by two simple examples. Finally we 
give a brief summary of the results. 



II. DIRAC PARTICLES IN ONE DIMENSION 

We work in (l+l)-dimensional space-time. The Dirac equation in an external vector field 
A v (t,x) reads 

(i-fD v -n)* = 0, (8) 

where \i is the mass of the particle, D v = d u + ieA u , e is the coupling constant, and summation 
over the repeated Greek index v {y — 0, 1) is implied. The 7 M are Dirac matrices satisfying 
the Clifford algebra: 

{y t ,7"} = 2<r, (9) 

where g^ u = diag(l, —1) is the Minkowskian metric. In this paper we only consider the zero 
component of A u , which is an even function of x, namely, we consider the special case where 

A 1 = 0, eA = V(x), (10) 

where V(x) is symmetric with respect to reflection: 
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V(-x) = V(x). 



(11) 



In this case we may set 

V(t,x) = e- iEt ip(x), 
and get a stationary equation for ip(x): 

Hip = Eip, 

where the Hamiltonian 

H = ap + 7°/ji + V(x), 



(12) 



(13) 



(14) 



where p = —id x is the momentun operator and a = 7°7 1 . 

To solve the Dirac equation (13) an explicit representation of the Dirac matrices is necessary. 
In one dimension this can be realized by the Pauli matrices: 
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7=0", 7 

In this representation a = —a 2 . We denote the two-component spinor ip(x) as 



(15) 



ip(x) = 



u[x) 



v(x) 



U\\X) 



u 2 (x) 



(16) 



where the second notation will be used only in a few occasions where summation over the 
spinor index is involved, then Eq. (13) can be explicitly written as a system of first-order 
differential equations for u and v : 



u' + (E + fi - V)v = 0, v' — (E — pi — V)u = 0. 



(17) 



Since V{x) is an even function of x, it is easy to show that if (u(x),v(x)) T is a solution of 
Eq. (17), then (u(—x), —v(—x)) T is also a solution with the same energy E, where r denotes 
matrix transposition. Thus the solutions of Eq. (17) can be chosen to have definite parities. 
Even-parity solutions are denoted by tp(x,+) and defined by the property of its components 
under reflection: 

u(-x,+) = u(x,+), v(—x, +) = —v(x, +). (18a) 
Odd-parity solutions are denoted by ip(x, — ) and defined by 



«(- 



-x. 



-) = -u{x, -), v(-x,-)=v{x,-). 



(18b) 



If (u(x), v (x)) T and (u(x), v(x)) T are both solutions to Eq. (17) with the same energy value 
E, it can be shown that uv — uv = constant. For bound states, all the functions u, v, u, and v 
must vanish at infinity, so the constant in the above equation is zero, and we have u/u = v/v. 
We denote this fraction by w, a function of x, and get u = wu, v = wv. Substituting this into 
Eq. (17), we have w'u = 0, w'v = 0. As a nontrivial solution, u and v cannot simultaneously 
vanish at any point, so we have w' = 0, or w = constant, which means that the two solutions 
are equivalent. Thus there is no degeneracy with bound states. 
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From the above discussions we conclude that a bound state solution of Eq. (17) must have 
definite parity. On the other hand, scattering states need not vanish at infinity. For a given 
energy, there are two linearly independent solutions. Both are physically acceptable. They do 
not necessarily have definite parities. However, they can be chosen such that one is of even 
parity and the other of odd. The above discussion is applicable to the case of free particles 
where V = and there is no bound state. 

For free particles, V = 0, then Eq. (17) can be easily solved. We have positive-energy 
solutions with E > /j, and negative-energy solutions with E < —fj,, both being scattering 
solutions. We define k = y/E 2 — /i 2 > 0, and denote positive- (negative-) energy solutions by 
the subscript k (—k), thus we have, say, 

E ±k = ±E k = iyW/i 2 . (19) 

The solutions will be given below. It is remarkable that when E = \i there is a nontrivial 
solution (not identically zero) with even-parity while when E = —/i there is one with odd- 
parity. These are uaually called half-bound states. 

Now we consider particles moving in the external symmetric potential V(x). In this section 
the potential need not be a cutoff one. We assume that V(x) decreases more rapidly than x~ 2 
when x — > oo, is less singular than x~ l when x — > 0, and is regular everywhere (except possibly 
at x — 0). It can be shown that E 2 > fi 2 gives scattering solutions while E 2 < fj 2 gives bound 
state solutions. (The critical energy states with E = ±pi will be discussed in the next section.) 
Scattering states will be denoted as above. Bound states with even (odd) parity will be denoted 
by a subscript k + («_) which takes discrete values. Scattering states are normalized by Dirac 
S functions. These are not necessary for the following discussions and are thus omitted. The 
orthonormal relation for bound states is given by 

I-+00 

(Vv ± ,-0k±) = / dxipl,(x)ip K± (x) = 5 K > K± , (20) 

J — oo ± 

which will be used below. Since the Hamiltonian (14) is an Hermitian operator, we have the 
completeness relation 

r dk ty k ( x , +)^ k (x', +) + -)v4(^> -) 

Jo 

+^_ fc (x, +)^_ k {x', +) + -)4>l k (x', -)] 

+ £ W {x)i>l + (x f ) + J2 1>*. (aOC = S(x - x>). (21) 

It should be remarked that all regular solutions should be included in the summation or 
integration. Thus half-bound states, when they appear, should be included, or the above 
equation does not hold. This can be verified straightforwardly in the free case where two 
half-bound states are present as pointed out above. The asymptotic forms of the various types 
of scattering sulutions are given by 



u± k (x, +) -> ±J cos[k\x\ + 77+(±£ fc )], (22a) 



v ±k (x,+) ^ J^±^e(x)sm[k\x\+r ]+ (±E k )}, (22b) 
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u± k (x,-) -> J^^e(x)sm[k\x\+r)_(±E k )}, (22c) 



U±fc (x,-) - Ty^^cos[A;|x| +77_(±f? fc )], (22d) 

when rr — > oo (+00 or —00), where e(x) = 1 (—1) for positive (negative) x, ?7 + (±i?fc) [r)-(±E k )] 
are the phase shifts of even-(odd-)parity solutions. The phase shifts depend on the sign as well 
as the magnitude of the energy, as Eq. (17) does. By setting all the phase shifts to zero, the 
right-hand side (rhs) of Eq. (22) gives exactly the solutions for free particles. Compared with 
the free solutions, the asymptotic forms in the external symmetric potential are distorted by 
the phase shifts. But it should be remarked that the normalization factors in Eq. (22) are the 
same as in the free case. 

It is well known that positive-(negative-)energy solutions correspond to particles (antipar- 
ticles) after second quantization. In the following we give some results for the scattering of 
positive-energy solutions by the symmetric potential V(x) described above. The results for 
the scattering of negative-energy solutions are similar. 

We denote the reflection amplitude by R k and the transmission one by T k . As a consequence 
of charge conservation, they satisfy 

|i4| 2 + |T fe | 2 = l. (23) 

Similar to the partial-wave method in three or two dimensions, R k and T k can be expressed in 
terms of the phase shifts. The result reads 

R k = l[ e *v + (E k ) _ e 2i V -(E k )^ = ie i [v+ (E k ) +v .(E k )] sin [ ?7+ ( jE ; fe ) _ v 4E k )], (24a) 
2 

T k = he***™ + e 2 *"-^*)] = e Hv + (E k )+v-(E k )] C0S [ v+ (E k ) - tj_ (£*)]■ (24b) 

These results obviously satisfy Eq. (23). They hold for both left-incident particles and right- 
incident ones since the potential is symmetric. If the potential is not symmetric, the reflection 
amplitude would depend on the direction of incidence. We can construct a scattering matrix 
and discuss various properties of it as in the nonrelativitic case [22]. However, we will not go 
further in this respect. In the following we will confine our discussions to symmetric potentials. 

From the above discussions, we see that information about the scattering process is con- 
tained in the phase shifts. The latter are determined by solving the system of equations (17) 
with the boundary conditions (22), and thus depend on the particular form of the potential. 
On the other hand, the total number of bound states with definite parities also depends on 
the particular form of the potential. The purpose of the Levinson theorem is to establish a 
relation between the phase shifts and the total number of bound states with a given parity. 



III. PHASE SHIFTS NEAR THRESHOLD 

In this section we discuss the behavior of the phase shifts r)±(E k ) and r)±(—E k ) near k = 
0. This is important in the development of the Levinson theorem, and is also helpful in 
understanding the result. 
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From now on we consider cutoff potentials that satisfy V(x) = when \x\ > a > 0. Such 
potentials will be denoted by V a (x) in the following. Since we will deal with solutions with 
definite parities, we need only consider the region x > 0. The solutions in the region x > a 
will be indicated by a superscript >. Scattering solutions are given by those on the rhs of 
Eq. (22) since V a (x) = in this region. Bound state sulutions will be discussed below. In the 
region x < a, the solutions are not explicitly available. First we consider the behavior of the 
solutions near x = and set appropriate boundary conditions. As assumed in Sec. II, V a (x) 
is less singular than x~ x when x — > 0. If V a (x) is regular at x — 0, the appropriate boundary 
conditions for Eq. (17) are obviously 

«(0,+) = l, u(0,+) = 0, (25a) 

u(0,-) = 0, u(0,-) = l. (25b) 

If V(x) behaves like A /x 1 ^ p ^ 2 when x — > + where A is a constant and < p < 2, it can 
be shown that the above boundary conditions remain valid. These boundary conditions are 
applicable to both scattering solutions and bound state solutions, with positive or negative 
energy values. Let us consider scattering solutions of Eq. (17) with positive energy E k that 
satisfy the boundary conditions (25). Note that Eq. (17) depends on k 2 rather than k, and the 
boundary conditions (25) are independent of k. Thus these solutions should only depend on 
k 2 as well. In other words, they are even functions of k. We will denote them by f±(x, k 2 ) and 
g±{x,k 2 ), where the subscript + (— ) indicates even (odd) parity. Thus scattering solutions 
with positive energy E k are given in the region x < a by 

u<(x, ±) = A ± (k)f ± (x, k 2 ), v<(x, ±) = A ± (k)g ± (x, k 2 ), (26) 

where the superscript < indicates the region x < a, A±(k) are constants which in general 
depend on k such that the solutions in the two regions can be appropriately connected. Ob- 
viously, u k (x, ±) and v k (x, ±) should be continuous at x — a, so that the probability density 
and the probability current density are continuous at the point. The phase shifts r]±(E k ) are 
determined by this condition. The results are given by 

/MQ-tanj; 1 + /M0 ^ f07 , 

tanr)+{Ek) = TTM^re tanr? -™ = tan e- mo ' (7a) 

where £, = ka and (3±(£) are defined by 



V E k -fj,f±{a, k 2 ) 

The above results show that the behavior of i]±(Ek) is determined by that of /3±(£) and 
ultimately by that of f±(a, k 2 ) and g±(a, k 2 ). The general dependence of f±(x, k 2 ) and g±(x, k 2 ) 
on k may be very complicated since Eq. (17) depends on k in a rather complicated way. 
Fortunately, only the property of f±(x,k 2 ) and g±(x,k 2 ) near k = is necessary for our 
purpose. 

We consider the limit E = E k — > /i (k — > 0) of Eq. (17). In this limit it takes the following 
form to the first order in k 2 : 



u (x 



+ P(x,k 2 )v(x) = 0, v'{x) + Q(x, k 2 )u{x) = 0, (28a) 



8 



where 

P(x, k 2 ) =2im+~V(x), Q(x,k 2 ) = -^- + V(x). (28b) 

The equations depend only on k 2 . Their solutions that satisfy the boundary conditions (25) 
should depend only on k 2 as well, since the boundary conditions do not depend on k. These 
solutions will be denoted by f±(x, k 2 ) and g±(x, k 2 ). Note that both P(x, k 2 ) and Q(x, k 2 ) are 
integral functions of k. Thus a theorem of Poincare tells us that f±(x, k 2 ) and g±(x, k 2 ), which 
satisfy k independent boundary conditions, are also integral functions of k. On the other hand, 
Eq. (17) coincides with Eq. (28) in the limit E = Ej, — > /i (k — > 0). Therefore, f±{x,k 2 ) 
and g±(x, k 2 ) must coincide with f±(x, k 2 ) and g±(x, k 2 ) respectively in the limit k — > 0, since 
they satisfy the same boundary conditions. Hence we conclude that f±(x,k 2 ) and g±(x,k 2 ) 
are analytic functions of k in the neighbourhood of k — 0. The above analysis holds regardless 
of whether the potential is cutoff or not, and the functions f±(x, k 2 ) and g±(x, k 2 ) are analytic 
in k near k = for any fixed x in the interval [0, +oo), not only in [0, a]. 

We have shown that f±(a, k 2 ) and g±(a, k 2 ) are even functions of k and analytic near k — 0. 
Therefore, when k — > 0, the leading term(s) for g±(a, k 2 )/ f±(a, k 2 ) must be given by one of 
the following forms: 

*tr 2l K ai + ateK (C-0) (29) 

where if, if, and if are natural numbers, of, af, af and af are nonzero constants. Conse- 
quently, the leading term for /3±(£) is given by one of the following forms: 

p±(0 - oc ± e i± -\ « ± r (2r± - 1} , (e-o), m 

where l± and l± are natural numbers, a± and a± are nonzero constants. Substituting this into 
Eq. (27a) we have 

tan^^)^^^ 1 or 6 + r (2p " +_1) , - 0), (31a) 

ta n?7 _(£; fe ) ^ &_r (2p " _1) or L^-\ (e-0), (31b) 

where p± and p± are natural numbers, b± and 6± are nonzero constants. Thus tanr/±(i?fc) — > 
or oo in the limit k — > or ^ — > 0, and r)±(fj,)/7i take integers or half integers. 

In order to study the threshold behaviour of r]±(Ek) more specifically, we consider bound 
state solutions of Eq. (17) with positive energy 

E = E X = yV 2 - A 2 , < A < /i. (32) 

The solutions satisfying the boundary conditions (25) in the region x < a will be denoted by 
F±(x, A 2 ) and G±(x, A 2 ). They are even functions of A as implied by the notations, since the 
equations are invariant under the change A — > —A and the boundary conditions are independent 
of A. Note that when —A 2 is replaced by k 2 , E x becomes E k and the solutions become the 
above scattering ones since the boundary conditions are the same. Thus we have the useful 
relations 

F±(x, -k 2 ) = f ± (x, k 2 ), G ± (x, -k 2 ) = g±(x, k 2 ). (33) 
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In the region x > a, the solutions are explicitly available since V a (x) = 0. They can be 
continuously connected to the interior solutions (for x < a) only for some specific values of 
A, which determine the discrete energy eigenvalues of bound states. This is quite different 
from the case of scattering states, where interior and exterior solutions can be continuously 
connected for any k if the phase shifts are chosen according to Eq. (27). On the other hand, 
given an energy value E\, a bound state with this energy eigenvalue exists only when the 
potential V a (x) has a specific form such that the solutions in the two regions can be connected 
continuously. Here we are interested in the case of the critical energy E — /j, (A = 0). The 
exterior solutions are then given by (up to a normalization factor) 

u>(x,±) = l, v>(x,±) = 0. (34) 

They can be continuously connected to the interior ones only when V a (x) takes some specific 
form such that 

G±(a,0)=0. (35) 

On account of Eq. (33), we conclude that critical energy solutions with E = jj, exist if and 
only if 

0±M) = O. (36) 

From Eq. (34) we see that the critical energy solutions, if exist, are not bound states since 
they are not normalizable. Nevertheless, the argument below Eq. (18b) applies since — 0. 
Thus the solution with E — /j, is not degenerate. In other words, the two solutions with E = \x 
and different parities cannot appear simultaneously for a given potential. As in the free case 
these solutions are called half-bound states. 

Now we easily realize that the first limit in Eq. (29) corresponds to the existence of critical 
energy states with E = \x. It can then be verified that this corresponds to the first case in Eqs. 
(30) and (31). The other cases in these equations correspond to the case without the above 
critical energy states. Thus we conclude that r] + (fi)/n [?7_ (//) / 7r] takes integers (half integers) 
when these exists a critical state with E = \x and even (odd) parity, otherwise it takes half 
integers (integers). It is easy to check that |T | = 1 when there exists a half-bound state with 
E = jj, (even or odd), otherwise T = [cf. Eq. (24)]. 

In the above we have analysed the threshold behaviour of rj±(E^) in detail. The threshold 
behaviour of r)±(—Ek) can be discussed in a parallel way. The results are given by 

t a n V+ (-E k ) - d + r {2q+ - 1] or d + e~ q+ ~\ (£-0), (37a) 

tanr7_(-£ fe ) -> d_f 2 *— 1 or tLr (2 *- _1) , (f -> 0), (37b) 

where q± and q± are natural numbers, d± and d± are nonzero constants. The first limit in 
Eq. (37a) or (37b) corresponds to the case when there exists a critical state with E = —\i and 
corresponding parity. Thus t) + (—(/,)/tt takes half integers (integers) when there 

exists a critical state with energy E = —[/, and even (odd) parity, otherwise it takes integers 
(half integers). 

The critical energy states with E = —/i, if exist, are given by (up to a normalization factor) 

v>Jx,±) = 0, v>(x,±) = l (38) 
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in the region x > a. These are also half-bound states. For a given potential, the two solutions 
with E = —fj, cannot apperas simultaneously. 

Using the threshold behaviours obtained above it can be easily verified that when k = Eq. 
(22) reduces to the above critical solutions or trivial solutions according as the corresponding 
critical solutions exist or not. 

It may be beneficial to discuss some specific model and to verify the above threshold 
behaviour of the phase shifts. A typical and simple cutoff potential is the square well potential 
with depth Vq and width 2a. This model can be solved explicitly, though the results are far 
from simple. Here we are interested only in scattering states and critical states. The conditions 
for the existence of critical energy states can be worked out explicitly. For scattering states, one 
can find the closed form for ta.nr]±(E k ) and tan^±(— Ek), and discuss the limit k — > 0. Since 
the calculations are somewhat tedious but straightforward, we will not give the details here. 
We just point out that the general results and conclusions obtained above are all confirmed by 
the simple model at hand. 



IV. THE LEVINSON THEOREM 

With the above preparations, we are now ready to establish the Levinson theorem. The 
theorem is developed on the basis of the completeness relations (21) and the threshold behavior 
of the phase shifts given in Eqs. (31) and (37). In addition, the fundamental equation (17) 
will be employed in the development of the theorem. 

Equation (21) is a matrix equation of which the rhs is a diagonal matrix. When written in 
matrix elements, it gives four equations. For free particles, the last two terms on the lhs are 
absent. We write down the first diagonal equation. Replacing x' by — x' in this equation and 
using Eq. (18) we get another. Taking the sum and the difference of these two equations we 
have 

r°° I 

/ dk [u° k (x, ±KV, ±) + Afo ±)u%W, ±)] = o M x - x ') ± S ( x + x % ( 39 ) 
Jo 2 

where the superscript indicates free particles. Carrying out the same procedure to the second 
diagonal equation we have 

roo I 

j o dk [v° k (x, ±K°V, ±) + v°_ k (x, ±)v°_* k (x', ±)] = -[S(x - x') =f 5(x + x% (40) 

The sum of Eqs. (39) and (40) gives 

r dk YlKs(x, ±)u° k * s (x', ±) + u°_ ks (x, ±)u°* ks (x', ±)] = 5(x - x f ), (41) 
Jo V 

where the second notation in Eq. (16) has been used and the spinor index is denoted by s. A 
similar result for the case with an external symmetric potential reads 

roo 

/ dk J2l u ks(x, ±)u* ks (x', ±) + u_ ks (x, ±)u*_ ks (x', ±)] 
Jo 

+EE^W< ±S M = ^-^')' ( 42 ) 

K± S 
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Now we subtract Eq. (42) from Eq. (41), then set x' = x, and integrate over x from — oo to 
+00, we arrive at 



n± = / dkM(±)^° k (±))-(M±),M±))} 

Jo 

+ / dk[(^ k (±),^_ k (±))-(^ k (±),^ k (±))}, (43) 
Jo 

where we have used Eq. (20) to get J2 K ±(' l f J K±,' l f J K±) — J2 K ± 1 — n ± where n + is the number 
of bound states with even (odd) parity. The inner products are defined by integrals similar to 
that in Eq. (20). From ort honor mal relations we have, say, (t/j k (±), t/j k (±)) = 5(0) ±8 (2k), and 
similarly for the other inner products in Eq. (43). These are infinities and should be treated 
very carefully. In order to avoid the difficulty of infiniteness, we define 



dx^\,(x,±)^ k (x,±), (44) 

-TO 



and obtain (^(±), ip k (±)) in the limit k! k and r — > 00. The other inner products in Eq. 
(43) will be treated in the same way. Using Eq. (17) it can be shown that 



(lp±k',^±k)r = ± p —(v* ±k ,U± k - U* ±k ,V± k ) 



ro 

(45) 

-ro 



which holds for both even-parity and odd-parity solutions. Using the asymptotic forms (22), 
and taking the limit k' — > k, we find 

GM+),lM+))n, = 7 + l ^f^ ± 2^E k Sin[2fcr ° + 2 ^ ±E ^ ( 46a ) 

{M-),M-))ro = 7 + l dV ~f k Ek) ^^- k si ^ 2kr ° + ^-( ±S *)1- (46b) 

For free particles, the corresponding results are obtained by setting the phase shifts to zero in 
the above equations. Obviously, the infiniteness in these results lies in the first term r /n in 
each equation when the limit ro — > 00 is taken. This disappears when we take the difference 
of Eq. (46) and the corresponding results for free particles. Using the well-known formulas 

hm = *(*), 

and g(k)5(k) = g(0)5(k) for any continuous function g(k), we obtain 



(^ ±k (+)^° ±k (+))r -(^±k(+)^± k (+)) 



ro 



l d V+ (±E k ) ± sin 2 v+{±fl)6{k) _p cos 2kr sin 2r) + (±E k ) , (47a) 

7r dk 2nkE k 

(r± k (-Ul k (-))r -(i>± k (-),i>± k (-)) 



l dr)-(±E k ) ^ sin 2 v {±fim) ± ^_ rCOs2 kr oS m2 V _(±E k ). (47b) 
7r dk 2irkE k 
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We regroup the four equations contained in Eq. (47) according to the sign of the energy instead 
of the parity, then integrate each group over k from to +00, and take the limit r — > 00, we 
have 

roo 



= -h±(/i)-^±(+oo)]±-sin 2 r ?± (/i)T7r lim / dk cos2kr , (48a) 



= -M-P) -ite(-oo)] =F isin 2 ^(-/i) ± f- lim fdfc 8 " 12 ^ ^ fc) cos2At , (48b) 

where we have used the integral j^ 00 dk5(k) = (1/2) dkS(k) = 1/2 since the Dirac 5 
function is an even function. So far in this section we have not cut off the potential. In the 
following we set V(x) = V a (x). Then we have the threshold behavior (31) and (37). The last 
term in Eq. (48a) can be decomposed into two integrals, the first from to e = + while 
the second from e to +00. The second integral vanishes in the limit r — > 00 since the factor 
cos2/cr oscillates very rapidly and the other factors in the integrand are finite. For the first 
integral, we have 

• a sm y m 2 tr , = i f dk s ^M = I r-^M) {m) 

kE k n Jo k n Jo 4 

since k < e is very small. On account of Eq. (31), we have sin 2r)±(E k ) — > c±!; 2r± ~ 1 (£ — > 0), 
where r± are natural numbers and c± are nonzero constants. Substituting into the above 
equation we have 



L 



;>!E^.£^il ^o, ( ^ 0+ ). (50) 



Thus the first integral vanishes as well. Then the last term in Eq. (48a) vanishes. By using Eq. 
(37) it can be shown in a similar way that the last term in Eq. (48b) also vanish. Substituting 
the results (48a, b), each without the last term, into Eq. (43) we arrive at 

[V±(v) ~ V±(+°°)} + [v±(-») - V±(-oo)} ± -[sin 2 r] ± (fi) - sin 2 rj±(-n)] = n ± n. (51) 

This is the Levinson theorem for Dirac particles in an external symmetric potential V a (x) in 
one dimension. It relates the phase shifts to the total number of bound states for each parity. 
In the next section we discuss some relevant problems and study two examples. Finally we 
summarize the results briefly. 



V. DISCUSSIONS 

1. Phase shifts at infinite momentum. It should be pointed out that there is no modulo-7r 
ambiguity in Eq. (51), because only difference of phase shifts at different momentums and 
trigonometric functions of the phase shifts are involved. However, it may be beneficial to 
appropriately define the phase shifts such that they can be determined uniquely. We consider 
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the potential 8V(x) where < 9 < 1 is a parameter independent of x, and V(x) is symmetric 
but not necessarily be V a (x). Positive-energy scattering solutions in this potential will be 
denoted by ip k {x, ±,#) and negative-energy ones by i[)- k (x, ±, 9). The corresponding phase 
shifts will be denoted by r]±(E k ,9) and r]±(—E k ,9) respectively. It is natural to define 

V± (E k ,0) = 0, V± (-E k ,0)=0, (52) 

since there is no potential in this case. We also require that rj±(E k ,9) and i]±(—E k ,9) 
be continuous functions of 9 for any finite k. Then the phase shifts r]±(E k ) = rj±(E k ,l), 
i]±(—E k ) = i]±(—E k ,l) in the potential V(x) are definitely defined. It should be remarked 
that the phase shifts at threshold are not continuous in 9, however (this was discussed in some 
detail in Ref. [18].). 

To determine the phase shifts at infinite momentum, we use Eq. (17) and Eq. (22) for two 
potentials 9V(x) and 9V(x). It can be shown that 

sm[A v+ (±E k , 9)] = T A9^(i> ±k (+, 9), Vi[> ±k (+, 9)), (53) 

where A9 = 9-9 and Ar] + (±E k , 9) = r}+(±E k , 9)-rj + (±E k , 9). In the limit -> 9 or A9 -> 0, 
the ip± k {x, +, 9) in Eq. (53) can be replaced by ip± k (x, +, 9), and the sine on the lhs can be 
replaced by its argument since the phase shifts are continuous in 9 as required above, so we 
have 

d v 4±J t ,e) = ^ (M+ , e) , VM+ , 0)) . ,54) 

Integrating over 9 from to 1 and using Eq. (52) we have 

V+ (±E k ) = =F— / d9 foM+, 9), V4> ±k (+, 9)). (55) 
k Jo 

For i]_(±E k ) we have a similar result. These results are practically not useful since the solutions 
on the rhs are not explicitly available in general. However, when k — > oo, we can ignore the 
potential 9V(x) in Eq. (17) since V(x) is not very singular at x = (less singular than x^ 1 as 
assumed) and is regular elsewhere. Then the solutions vp± k (x, +, 9) in Eq. (55) can be replaced 
by the free ones i[)± k (x, +) and we have 

V+{±oo) = =f / dxV(x). (56) 
J o 

The result for r]_(±oo) is the same. The integral in the above equation converges because 
V(x) decreases more rapidly than x~ 2 when x — > oo and is less singular than x^ 1 when x — > 0, 
as assumed in Sec. II. These results are similar to those obtained in three dimensions [9, 16, 
17] and two dimensions [18]. Of course, they hold in the special case V(x) = V a (x). As a 
consequence, we have 

^±(+oo) +v±{-oo) = 0, (57) 

and the Levinson theorem (51) reduces to the form of Eq. (7). 

We have calculated i] + (±oo) and ^_(±oo) exactly for the square well potential, and the 
results (56) are confirmed. In the following we will see that Eq. (56) does not hold for the 5 
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potential well, however. This is because the 5 potential well is more singular than the type we 
have assumed. Indeed, if V(x) is less singular than x^ 1 at x — 0, or behaves like Aqx~ 1+p ^ 2 , 
we have J £ dx V(x) = 2A e p ^ 2 /p — > (e — > + ). If V(x) = — £/o<5(x), however, we have 
Jg diy(x) = —Uq/2. Thus the 5 potential well is more singular. It is expected that Eq. (56) 
remains correct as long as V(x) belongs to the type we assumed, i.e., less singular than x~ x at 
x = 0. 

2. Verification of the theorem. To examine the Levinson theorem one should choose a 
simple potential such that n± and all phase shifts at zero momentum can be worked out 
explicitly. This is not available even for the square well potential. Although closed forms for 
tanr)±(Ek) and tan?7±(— Ek) can be obtained, they depend on k in a very complicated way. 
To determine r)±(ii) and r)±(—ii), numerical calculations are necessary. The transcendental 
equations for the energy levels of bound states are also complicated. Perhaps the simplest 
potential is the 5 potential well V(x) = —UqS(x), where Uq > is a dimensionless parameter. 
We have pointed out above that this potential does not belong to the type we have assumed 
in developing the Levinson theorem. However, threshold behavior of the phase shifts similar 
to Eqs. (31) and (37) can be explicitly shown. So the Levinson theorem should remain correct 
in this case. As the calculations are simple, we only give the results. We have n + — 1, n_ = 0. 
The phase shifts at infinite momentum are different from those given by Eq. (56): 

r] + (±oo) = iarctan-^, (58) 

where arctan(?7 /2) G (0, tc/2) is the principal value. A similar result holds for ?7_(±oo). Thus 
Eq. (57) remains valid in this case though Eq. (56) does not, and the reduced form of the 
Levinson theorem (7) is expected to hold. In fact, the phase shifts at threshold can be found 
to be 

71 71 

V+(V>) = 2> ^("^ = °> ^-(^) = °> = - g' ( 59 ) 

and it is easy to verify that Eq. (7) is satisfied. 

For a S potential barrier V(x) = Uo5(x) where Uq > 0, it can be shown that n + = 0, 
n_ = 1, rf + (fx) = — 7r/2, r] + (—fi) = 0, rf-(fx) = 0, rf-(—fx) = 7r/2, and Eq. (7) is satisfied as 
well. 

A less simple example is the double 5 potential wells V(x) = —Uq[5(x — a) + 8(x + a)], 
where Uq > is dimensionless. This example is somewhat more substantial since half-bound 
states may be involved, and both n± and the phase shifts depend on the value of Uq. Since 
the calculations are not difficult but lengthy, we will not give the details. We only point out 
that in this case the Levinson theorem is confirmed once again. 

3. Summary. In this paper we study Dirac particles in one-dimensional symmetric po- 
tentials that decrease more rapidly than x~ 2 when x — > oo and are less singular than x~ x 
when x — > 0. The properties of bound state and scattering solutions are discussed. For cutoff 
potentials the threshold behaviours of the phase shifts are studied in detail, and the Levinson 
theorem (7) is established, which connects the phase shifts with the total number of bound 
states for each parity. Two simple examples are discussed and the Levinson theorem is ver- 
ified explicitly. A mathematically rigorous extension of the Levinson theorem to non-cutoff 
potentials is still not available. However, one may expect that the theorem remains correct 
for short-range non-cutoff potentials where the asymptotic forms (24) for scattering solutions 
hold and the total number of bound states is finite. 
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